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Abstract 

Thermal two-loop QCD corrections associated with light stops have a dramatic ef- 
fect on the strength of the MSSM electroweak phase transition, making it more strongly 
first order as required for the viability of electroweak baryogenesis. We perform a per- 
turbative analysis of the transition strength in this model, including these important 
contributions, extending previous work to arbitrary values of the pseudoscalar Higgs 
boson mass, uia- We find a strong enough transition in a region with 2 ^ tan/? ^ 4 
and niA ~ 120 GeV, a light Higgs boson with nearly standard couplings, and mass 
below 85 GeV within the reach of LEP II, and one stop not much heavier than the top 
quark. In addition, we give a qualitative discussion of the parameter space dependence 
of the transition strength and comment on the possibility that the transition turns to 
a crossover for sufficiently large Higgs masses. 
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1 Introduction 



It is by now well established that successful electroweak baryogenesis calls for new physics 
at the Fermi scale. In the minimal Standard Model (SM) it is not possible to generate a 
sufficient baryon asymmetry at the electroweak phase transition (EWPHT) because the 
necessary CP violation is far too small. In addition, the phase transition is at best weakly 
first order for realistic values of the Higgs mass, and any net B + L number created would be 
subsequently erased by unsuppressed sphaleron processes in the broken phase. This failure 
could be an indication that some type oi B — L violating new physics is at work above 
the electroweak scale (see e.g. [^). Although that is an interesting possibility, we are not 
really forced to give up the beautiful idea of electroweak baryogenesis, because physics at 
the 100 GeV scale may not be described by the SM. Supersymmetry (SUSY) 0, the most 
promising candidate for physics beyond the SM, is expected to be relevant at such energies. 
Thus, it is natural to examine the prospects for electroweak baryogenesis in models of low- 
energy SUSY. This is particularly interesting for two reasons; first, because low-energy SUSY 
is well motivated for reasons not related to the matter asymmetry problem. Second, the 
parameter regions where electroweak baryogenesis may be successful are in the reach of the 
present generation of colliders so that the viability of the scenario would be tested in a not 
too distant future. 

Supersymmetric electroweak baryogenesis has been the subject of many studies in recent 
years. We restrict ourselves to the simplest realization of low-energy SUSY, the Minimal 
Supersymmetric Standard Model (MSSM). Non-minimal modelsQ have also been considered 

i- 

In SUSY models the presence of extra sources of CP violation besides the Kobayashi- 
Maskawa angle in the SM can be enough for electroweak baryogenesis and imply CP violating 
signals on the verge of being seen in upcoming experiments. The observed baryon asymmetry 
can be accounted for if the new CP violating phases are greater than lO^^^^^^ and some 
superpartners are light enough^ [0 (see also p, |, |^). In addition, the presence of two 
Higgs doublets makes possible the spontaneous violation of CP in the Higgs sector at finite 
temperature |T^ (spontaneous breaking does not occur at T = for realistic values of 
the parameters of the model This mechanism takes place for small values of the 

pseudoscalar Higgs boson mass uia and large values of tan/? (the ratio of Higgs vevs), and 
could play an important role for electroweak baryogenesis. 

The requirement of a sufficiently strong first-order transition to avoid B + L sphaleron 
erasure of the asymmetry sets a strong upper limit on the mass of some Higgs boson in the 
theory. Therefore, it is crucial to determine the strength of the electroweak transition and 
identify the regions of parameter space where it can be strong enough. More precisely, the 

^In models without i?-parity conservation, baryogenesis before the EWPHT faces the problem of baryon 
depletion by the combined action of sphalerons and B ~ L violating processes. Electroweak baryogenesis 
offers a particularly appealing solution to this problem, although other possible solutions have been proposed 

i- 

^These results are affected if the rate of baryon number violation in the high temperature symmetric 
phase is smaller than it was previously thought to be (see [jll], |l^ p^). 
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requirement is |^ 



Tr 



> 1, 



(1) 



where Tc is the critical temperature of the transition, defined by the coexistence of two de- 
and V = 



generate mmima, 



+ V2, the order parameter of the transition, is the vacuum 



expectation value (vev) of the Higgs field driving electroweak symmetry breaking (normal- 
ized to f = 246 GeV at T = 0). The first studies of this problem in the MSSM have already 
pointed out that (|1|) would require light stops in the spectrum (that is, with masses compa- 
rable to the critical temperature of the transition). In ref. it seemed inplausible that 
these constraints were satisfied, so MSSM electroweak baryogenesis was considered unlikely. 
On the other hand, ref. [0 found a sizeable region of parameter space with light stops 
and large mass for the pseudoscalar Higgs boson, in which the transition was strong 
enough for baryogenesis. However, a careful treatment of the LEP bounds on the Higgs 
mass for large tua reduces that region significantly. Both fl^ and ||18| studied the one-loop 
effective potential perturbatively at finite temperature without resummation of higher loop 
contributions. 

A more careful analysis of this potential, including resummation of the so-called Daisy 
diagrams which have an important effect, was carried out in [jl9[ (for large m^) and ||2^ (for 
arbitrary m^). A small region of parameter space in which the transition is strong enough 
remained after imposing experimental constraints. This region corresponded to large rriA 
(which in turn determines a Higgs spectrum with only a light Higgs boson), small tan/? 
(tan/3 ~ 2), light stops compatible with experimental limits and with negligible L-R mixing, 
heavy gluinos, and the light Higgs barely above the LEP lower limit {rrih > 65 GeV at that 
time). The rest of the supersymmetric particles do not influence much the transition so 
that their spectrum is not constrained, although the general tendency is to prefer heavy 
superpartners. Even if this is a signiflcant improvement with respect to the SM, the window 
for baryogenesis was clearly small and is already excluded after the improvement in the Higgs 
mass bound from LEP {rrih > 70.7 GeV from ALEPH 0]). 

22, 23 1 . In ref. |2^, the region of light masses for 



The subject was revived recently by refs. 



tn was increased by taking negative values of its soft-mass squared, mfj < 0. Negative values 
of mfj are associated with the existence of dangerous colour-breaking minima in the MSSM 



multi-scalar potential. While in |]T9|, |2^ the lower limit on was set to zero, ||22[ shows 



that moderately negative values of mfj can still be cosmologically allowed; the electroweak 
minimum is the deepest one, and is the one populated when the Universe cools down. By 
moving into a region with smaller m^^ (while keeping negligible mixing), the electroweak 
transition is stronger and sphaleron erasure can be avoided (see also |2^), thus enlarging 
the region where electroweak baryogenesis could take place. 



In ref. 



it was shown that two-loop QCD corrections associated with light stops 
dramatically affect the quantitative behaviour of the effective potential. Inclusion of these 
contributions is then mandatory for a more precise study of the transition, and these correc- 
tions make it strong enough for baryogenesis even if mfj > 0. Therefore, this effect improves 
the situation in all the regions of parameter space with light stops. 

Both ||2^ and focused on the large m^-small tan /5 region of parameter space preferred 
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for baryogenesis and used perturbative calculations of the temperature dependent effective 
potential. As discussed in the perturbative treatment is expected to give reliable 

results for Higgs masses larger than in the SM case (where the limit for the validity of 
perturbative calculations was roughly around mw). A very useful complementary approach 



was followed in refs. where 3d effective theories |2^, |3y] were constructed to study 

the high temperature phase of the MSSM and, in particular, the nature of the electroweak 
transition. For not too small mfj, the relevant 3d effective theory is well approximated by 
an SU{2) + Higgs model (as in the SM but with different couplings) for which there exist 
non-perturbative lattice studies. The results of these analyses are in reasonable agreement]^ 
with those of |2^, A region with a strong transition is confirmed for light stops (in 

particular light tj?), small stop mixing, small tan/5, and Higgs masses up to about 80 GeV. 

In any case, it is important to keep in mind that for small values of mu (the interesting 
region for baryogenesis) the relevant 3d effective theory is no longer SU (2) + Higgs, but tji 
and gluons should also be taken into account. The construction of such an effective theory 
is straightforward and was sketched in |^ , but Monte Carlo simulations for that theory are 
not available. Another possible shortcoming of the existent 3d reduction studies may show 
up in the small rriA region where the two Higgs doublets are light. In such a case the effective 
theory may require the presence of both Higgses; it may not be justifiable to integrate out 
the heavier of the two. This point is discussed in more detail in the next sections. 

The purpose of this paper is to further continue the perturbative exploration of the MSSM 
baryogenesis window in the region of small and moderate m^, including the important two- 
loop corrections. As we will show, smaller ttia always corresponds to a weakening of the 
transition so that (for fixed values of other parameters) there is a lower bound on rriA below 
which the transition is too weak for baryogenesis. This lower limit decreases for smaller 
tan p. 

The small tua case is interesting for several reasons. First, it can be preferred in some 
mechanism for the generation of baryon number [0. Second, in the case of small m^-large 
tan (3, the spontaneous CP violation mechanism of ref. already mentioned could play an 
important role in baryogenesis. Unfortunately, we do not find any allowed region for large 
values of tan/3. Finally, the 3d approach in this region may be less straightforward than for 
large tjia- In particular, contact with the Monte Carlo lattice results for the SU{2)+ Higgs 
theory could be hampered by the necessity of keeping both Higgses in the effective theory. 

In this paper we are going to make a purely perturbative analysis of the transition, 
so we take some time to discuss in section 2 the range of parameter space where non- 
perturbative effects start to play a non-negligible role in the phase transition. In analogy 
with the SM case, we expect that this range is associated with the change in the nature of 
the EWPHT from first order to a crossover. We discuss in section 3 the expectations for 
a strong phase transition in the (m^,tan/3) parameter space making a simple qualitative 
analysis of the effective potential, both at zero temperature and including the dominant 
temperature corrections. In section 4, we write the high temperature effective potential at 



Some of the results obtained in p8| were artifacts of an implicit expansion in the stop mixing over the 
Debye masses and disappear when the convergence of this expansion is under control ||3l|] . In particular, the 
baryogenesis region of small m^-large tan/3 found in does not subsist. 
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one-loop order with Daisy resummation, concentrating on the small itla region. In section 5, 
we consider the effect of dominant contributions to the two-loop resummed potential. We 
present our results in section 6 and our conclusions in section 7. In appendix A, we collect 
field dependent masses and mixing angles (both at zero and finite T) for the particles relevant 
to our study. Appendix B gives the lengthy expression for the non-expanded full two- 
loop potential. Finally, appendix C contains the dominant contributions to the two-loop 
resummed potential in a high T expansion. 



2 Crossover in the MSSM 

In the MSSM with light stops, the expansion parameter of the resummed perturbation theory 
at finite temperature is eMssM ^ hf/\ rather than esM ~ 9^/^ HH HI HI) where ht is the 



top Yukawa coupling, g the SU{2)l gauge coupling and A the quartic Higgs boson coupling. 
Perturbation theory is expected to be applicable to the study of the EWPHT in the MSSM 
for a range of Higgs masses wider than in the SM, with a critical upper limit governed by mt 
instead of mw- We present a rough estimate of a related quantity, the critical Higgs mass 
beyond which the phase transition changes from first-order to a crossover. In the SM, it 
is possible to estimate this mass analytically by imposing the equality of the transverse W 
mass in the broken and symmetric phases at the critical temperature The argument 
arises naturally in the context of the 3d reduced theory but we use 4d quantities in what 
follows. The mass in the broken phase is given by 



mwiTc)i,^ = ^gv{T^), (2) 



where the vev v{Tc) can be obtained from the effective potential 

V ^ lm\T)v' + lx{T)v' - ^gV. (3) 

2 8 loTT 

The mass in the symmetric phase, the magnetic mass, is of non-perturbative origin and can 
be estimated by solving a set of gap equations [PSl ^M. It is 



mw{Tc),y^ = Cg'T,, (4) 

where C = 0.28 in the Standard Model (it is basically equal to its value in the SU (2)-(t 
model). Equating (|^) and (|^) and solving for A one arrives at the critical mass for the onset 
of crossover 

ml = ^l^rnw ~ 74 GeV . (5) 

This naive estimate is close to the numerical result of lattice simulations |3^ which give 
ml - 80 GeV. 
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In the MSSM case with large and hght stops, we add to the potential the stop 
contribution, roughly approximated by 



SV ^ -r^-j=hlsM^\ (6) 

where ht^sM = ht sin /5, and the constant r parameterizes the effective strength of the stop 
corrections. It is normalized in such a way that it would be equal to 1 if all screening effects 
from soft and thermal masses were negligible. In realistic cases r ^ 1. Heavy soft masses or 
light gluinos (which give a sizeable contribution to the stop thermal mass) would decrease r. 
On the other hand, a negative value for m\j tends to increase r. From the numerical study 
of the transition, we can estimate that r ~ 0.2 when stops affect the transition sizeably 
without conflicting with experimental constraints. 
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Figure 1: Critical Higgs boson mass for crossover in the MSSM [as estimated by formula (^]. The 
parameter r measures the effect of stops on the transition: r = corresponds to the pure SM limit and 
realistic values of r should not be larger than 0.2 — 0.3. 

The magnetic mass does not depend on the presence of stops (in the same way that it is 
not sensitive to the Higgs mass in the SM) and it should be well approximated by the same 
expression (^). The following estimate for in the MSSM follows 



ml = mw 



rrit 



AttC 



w. 



(7) 
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Eq. d^) is plotted as a function of the parameter r in figure 1. Due to the factor, 
grows very rapidly with r; for example, r ~ 0.1 corresponds to ml ~ 100 GeV. Below the 
curve the electroweak phase transition is first-order, and weakens as the curve is approached. 
Above the line there is not a phase transition, but an analytic crossover. 

To interpret figure 1 correctly one needs to take into account that a change in the masses 
of the stops affects r directly, and ruh through radiative corrections. Moreover, we have 
obtained eq. (0) for the case of large m^, when A in @ can be directly related to the 
light Higgs boson mass. For lower rriA one can still obtain a critical value A'^ for A, but its 
relation to m% involves the mixing angle between the two Higgs scalars (and < X'^v). 
Also, the important effect of two-loop corrections on the strength of the transition cannot 
be accounted for by a simple redefinition of r. The results presented in figure 1 are thus 
a rough estimate of the critical mass values expected. An alternative way to describe the 
emergence of crossover is to equate and (0) to obtain v/Tc ~ 2Cg ~ 0.4, and interpret 
this as the critical value for the jump in the order parameter as computed in perturbation 
theory. When perturbatively the transition is so weak, non-perturbative effects can no longer 
be neglected and they would eventually change the transition from first order to a crossover. 
This definition can be used now with the two-loop corrected v/Tc. As we are primarily 
motivated by the possibility of electroweak baryogenesis, which requires a strong transition, 
we do not attempt here a more precise determination of along these lines. 



3 Parameter Space 

Light stops increase the strength of the electroweak phase transition in the MSSM. Whether 
this effect is large enough to permit baryogenesis depends on the details of the Higgs sector 
at T = 0. At tree-level, this sector is completely determined by two parameters: the 
pseudoscalar mass ttia, and tan/3, the ratio of Higgs vevs at zero T. In this section we 
explore the plane (mA,tan/3) showing the expectations for a strong phase transition. 
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Large (m^ > mz,Tc). 



In this case only one Higgs doublet remains at the electroweak scale and is related to 
symmetry breaking, while the other is heavy and decouples from the problem. The Higgs 
potential relevant to the study of the electroweak transition is SM-like, V{ip), with = 
ipi cos P + (p2 sin P, but the light Higgs mass is not a free parameter as in the SM and it is 
determined at one-loop to be 



ml = ml cos' 2/5 + A !^ log + /(A, + /./ tan/5). (8) 

The logarithmic radiative correction increases with increasing stop masses and can be size- 
able. The last term in (H) depends on the stop mixing through the combination Xt = 
At + fi/ tan/3 and reads: 



3 m 



fiXt) = T-A^t [2h{ml,ml)+X'^f{mlmi)\ , (9) 



47]-2 y2 
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where 



and 



h{a, h) 



1 a 

^^°St, 

a — 



by 



2 - 



6 a 



(10) 



'111 



The correction (^) is zero if the stop mixing is neghgible ("zero mixing" case), reaches a 
positive maximum for large values of the mixing ("maximal mixing" case), then drops (even 
getting negative) if the mixing is further increased ("extreme mixing" case). 

As in the SM case, the phase transition becomes weaker for larger m/i, so the best region 
for baryogenesis always corresponds to a region with lightest possible Higgs massQ. From 
(H), we see that the region preferred to get a strong transition is one with small tan/5, and 
light stops with negligible mixing (the case of a light Higgs due to extreme mixing will be 
discussed later). Of course, it is not enough to have a light Higgs if no effect beyond the 
SM helps to strengthen the transition. There are two benefits of having light stops; first, 
the radiative corrections to the Higgs mass are smaller for lighter stops so that the Higgs 
boson can be lighter. Second, direct stop contributions to the T 7^ potential enhance the 
transition strength (for small mixing). These ingredients describe the well known region 
preferred for baryogenesis. We now examine other possible effects to see if we can find 
additional regions supporting a strong phase transition. 

• Intermediate mass stops 

If the stops are very heavy compared to the transition temperature, their thermal contri- 
butions are Boltzmann suppressed and the SM case for the corresponding rrih is recovered. 
On the other hand, moderately heavy stops can still influence the transition via a T = 
effect described in [^] (see also the discussion in fl^l)- If we neglect all couplings other than 
ht in the loop corrections, the one-loop T = potential can be written approximately as 



+ 



6 



647r2 



4 

rrir 

ti 



log 



ml 




2^^(logf-{ 



, (12) 



where m^^ ^ '^i field-dependent quantities and can be found in appendix A (making 
the replacement (fi — > (pcos(3, <^2 ~^ fsinf], appropriate in the large limit). When 
stops are heavy (for example, with equal soft masses rriQ = mu = Ms) compared to mj, it 
is possible to expand their contribution to the potential in powers of (p/Mg. If we express 
and A of eq. (|T2|) in terms of the radiatively corrected vev v and mass rrih, the potential 
takes the simple form 



V(v^'-^')^ + A\4(^)+A\^-(^), 
4 



(13) 



'*The fact that the light Higgs is SM-hke imphes that the usual LEP bound on its mass {irih > 70.7 GeV) 
also applies; for large tjia there is no possibility of accessing lower values of mh to strengthen the transition. 
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with 

4 2 

^Vt{^) = 

and 

4 2/ 2 4 




(14) 



^^'^^^ = r ~ V + + ^(^7M|). (15) 



Here 



2M|' 



(16) 



^^t,ti.v) satisfy AK' = AK" = at = f so that they do not affect v or m/^. Although /^Vi 
shows the decouphng exphcitly, it represents a deviation from the pure SM potential for not 
so large values of Ms- This effect corrects the transition strength for a fixed Higgs mass. For 
example, if Ms = 500 GeV the increase in v{Tc)/Tc with respect to the SM with the same 
Higgs mass can be estimated to be a factor 1 + 0.1k larger. This is a 10% effect if the stop 
mixing is negligible. Unfortunately, this modest increase for heavy stops is counterbalanced 
by the direct increase of through radiative corrections so there is no net gain, and larger 
stop masses always give a weaker transition. On the other hand, it is apparent from ([T6|) 
that stop mixing tends to weaken the enhancement of v{Tc)/Tc, eventually having a negative 
effectg. 

• Extreme stop mixing 

In this case, one may hope to get small masses for the Higgs boson via a negative radiative 
correction from the last term in To reach this region either one or both stops should be 
heavy, otherwise the large mixing would drive the lighter stop mass to an imaginary value. 
When both stops are heavy, the situation is described by the preceding paragraph. If one 
stop remains light (say in, to avoid problems with Ap), its mass squared takes the form 



K = < + nh't sin^ P^'{l-^]+ Dni^), (17) 



1 / 

in which Dptif) is the D term contribution, and niu {mq) is the tn (ti) soft mass. Eq. (p!7| ) 
shows that the coupling to the Higgs field ip is reduced by the mixing, and the influence of 
this stop on the transition diminishes [^. Thus, no improvement is expected over the SM 
in this region. 

• Stop mixing two-loop effects 

These corrections, calculated in [^, can have two different effects. First, the stop mixing 



angle at appears in the stop contributions. This changes the corrections slightly but no 
dramatic enhancement appears. The second effect, potentially more important, arises from 



^The rest of squarks and sleptons give only a small correction (see ||l7|]) through this effect. On the other 
hand, if they are lighter and in thermal equilibrium they raise the screening masses of stops and thus have 
a negative effect on the transition strength (see however p^). 
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a set of corrections that depend on the new trihnear couphng h ~ t — t. The dominant 
contribution of setting sun diagrams involving this couphng is ||23|| : 



N 



og 2^ + og 



^ + log ^ 

Barred scalar masses in this expression include effects from thermal screening. We have 
omitted in (|18|) other contributions proportional to gauge couplings or /ife, as well as terms 
of the type \og[{mh + 2mfJ/(m/i + 2mi,y\, that have a small effect on v/T^. (a more refined 



approximation can be obtained from ref. 



It is a simple exercise to show that a contribution to the potential of the form 5V = 
log[((y9^ + n)/n] increases v/Tc for positive K (see ||38|). The contribution (0), which can 



be approximated by an expression of this form, tends then to raise v/Tc- However, keeping in 
mind that Xt cannot be made much larger than the scale T, rough numerical estimates show 
that this two-loop enhancement of v/T^. is below the few percent level, while the negative 
effect of a non zero Xt at one- loop level [O] is quite significant. 



2) Small niA {jtia ~ mz,Tc). 

For low rriA the transition is not forced to proceed along the fixed direction (p2/'^i = tan 
Consider the (T ^ 0) mass matrix at the origin ipi = = and define the field-direction 
of the lowest eigenvalue as = ipiCos6{T) + ip2sm6{T). The behaviour of 0{T) has been 
studied in refs. 0. At T = it is given by 

tan 2/3 , , 

m^ + M| + AmVcos2/3 ^ ^ 

where Am^ > is a one-loop correction term that can be found in We see that 

tan^ > tan/3 and for sufficiently low rriA, ta.n9 ^ tan/3. When T increases, the radiative 
corrections in (|19|) receive a negative contribution that grows like h'^T'^, which decreases 
tan^(T). The critical temperature for the transition is usually reached for tan ^(Tc) 3> tan/3 
(see 

In the low m^-low tan/3 (~ 1) region, 9(Tc) could be below tan/3 In general, 

tan^(Tc) is approximately equal to ip2(Tc) / ipi(Tc) with the larger differences being ex- 
pected for low tan/3. If tan ^(Tc) and (p2(Tc) / (pi(Tc) differ significantly, the usual procedure 
of integrating out the heavy Higgs at Tc may not be justified in the construction of 3d ef- 
fective theories, at least to study v/Tc. However, as already noted in this discrepancy 
is most significant for values of rriA and tan/3 excluded experimentally. Nevertheless, it is 
important to keep this possible complication in mind when constructing 3d effective theories 
for models with a non-minimal Higgs sector. 



^In the language of 3d effective theories, the heavy Higgs field that can be integrated out at the moment 
of the transition is basically Hi for small m^. 
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In the most common case, 0{Tc) ~ (p2{Tc) / ^pi{Tc) , and we can safely assume that the 
transition takes place along that direction. In particular, we note that the Higgs quartic 
coupling relevant to determine the strength of the transition is the quartic self coupling of 
that particular direction. In other words, the Higgs mass whose magnitude controls the 
transition strength is some effective mass for excitations along the field direction excited 
in the transition^. In principle, this is not the mass of any of the physical Higgses. It is 
straightforward to obtain an expression for this effective mass which is the one written in 
eq. (H) with the replacement f3 0(Tc). At low rriA, the physical mass of the lightest Higgs is 
below (^) and the corresponding eigenstate is associated with excitations in a field direction 
different than 6{Tc). 

If one starts with some v/Tc for large and fixed tan/5, the two effects just described 
cooperate to decrease v/T^ as rriA lowers; lower values of rriA increase tan 6'(Tc), which probes 
larger effective masses and thus weakens the transition. The maximal value of (§) is already 
saturated for large tan/3, so tan 6'(Tc) ^ tan/? does not make much difference in that case: 
in the large tan/5 regime the transition tends to have v/T^ small and independent of tan/?. 



This feature is actually observed in refs. |26, 27|. The only low niA region where one can 



hope for a strong transition is once again the region with low tan/?, although it is expected 
to be reduced with respect to the large ttia case. 

The qualitative behaviour described in the previous discussions is confirmed by numerical 
computations of v/Tc. In the rest of the paper, we focus on the small region and compute 
v/Tc by studying the effective potential up to two loop order in resummed perturbation 
theory. 



4 One-loop resummed potential 

In this paper, we assume that the only particles present at the energy scale of the electroweak 
phase transition are the SM particles plus the extra MSSM Higgs doublet and squarks of the 
third generation [Ql = itL^bi) and tij,&_R]. The masses of the rest of the supersymmetric 
particles are assumed to be large compared to the critical temperature of the transition 
and their contributions are Boltzmann-suppressed. This is just a simplifying assumption 
and our results are qualitatively the same for more general spectra. The most important 
point is that gluinos should be heavy and Boltzmann-decoupled. If not, they would give 
significant contributions to the thermal masses for stops [T^, , making them heavier and 
weakening the phase transition. The details of the transition depend weakly on the rest of 
the spectrum. Our model gives a fair estimate of the phase transition strength for realistic 
supersymmetric spectra (in particular models with a neutral LSP), with the only condition 
that gluinos should be heavy. 

For general values of m^, the effective potential is a function of two Higgs background 
fields: = {H'()\/2 and = {H2)y/2, in which are the neutral components of the two 
MSSM Higgs doublets. At T = 0, ip2/^i = V2/V1 = tan/? and + = (246 GeVf. 

similar observation is helpful to estimate rather accurately the sphaleron mass in the MSSM [p9[. 
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The tree-level effective potential is 



1 1 g"^ + g''^ ■ 



2\2 



(20) 



with the tree-level quartic couplings fixed by SUSY in terms of gauge couplings as shown. 

At one-loop (with resummation of Daisy diagrams), the effective potential (written in 
MS scheme and 't Hooft-Landau gauge) receives the contribution: 



Vl{ipi,ip2) 



1 

2^ 



mi{^i,^2) 



(21) 



The field-dependent masses m'^{ipi,ip2) for gauge bosons, top and bottom, Higgses and 
third generation of squarks are given in appendix A, both at zero temperature and in the 
thermal plasma. The n^'s are the number of degrees of freedom and can also be found in the 
appendix A, Q is the renormalization scale, and Ci = 5/6 (3/2) for gauge bosons (scalars 
and fermions). The mass parameters 2 can be traded by by minimizing the T = 
one-loop potential |2y]. The Jj's give the finite temperature effects; the tilde denotes that 
daisy resummation has been performed in the Jj's. Here Jj = J± is the free energy of an 
ideal gas of particles of mass mi{(pi, ip2), 



Mmyr^) = dxx' log (^1 T e-v/^^+^^) 



f22) 



where +(— ) is for bosons (fermions). No Daisy resummation is needed for fermions and 
J_ = J_ . In the scalar sector, we choose to perform this resummation not only on the zero 
Matsubara modes but on all of them Bl|, so that 



(scalar)/ 2\ 



J4 



(mj) 



(23) 



where the rnf are the thermally corrected masses as given in appendix A. We perform daisy 
resummation for gauge bosons by screening the n = modes of longitudinal components 
Wl, Zl, -fL, gL- 



i,long) 



J4 



71 m, 
6 



3 

iL 



T3 



(24) 



where rriiL are the Debye masses (to be found in appendix A). In eq. (|2^ the n = (cubic) 
contribution to J+ is subtracted out and replaced by a similar term with the thermally 
corrected mass. Transverse modes are not screened at leading order. 

This potential was studied in [^, where it was shown that the presence of light stops 
can influence the strength of the electroweak phase transition through their contribution to 
the cubic n = terms. The final region in parameter space where v/Tc > 1 is determined 
by the interplay between two opposite effects; soft and thermal screening masses tend to 
decrease v/T^ by screening the pure cubic behaviour of m?, while a large Yukawa coupling 
tends to increase v/Tc. The numerical study of the potential shows that in some region of 
parameter space the cubic term from stops can dominate the electroweak phase transition. 
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However, the only region where this transition was strong enough for baryogenesis was 
hmited to the large ttia, small tan/3, negligible stop mixing regime. For example, with a 
top quark pole mass Mf = 156 Gel^J^ stop soft masses mq = 70 GeV, rriij = 0, zero stop 
mixing At = fi = and tan/5 = 2.5 the ratio v/T^, which is ~ 1 for large tua, drops below 



0.5 for rriA = 50 GeV. This is shown in curve (6) of figure 2. Here v = yvf + V2, and is 
defined by the coexistence of two degenerate minima in the T-dependent effective potential, 
one at {{pi,ip2) = (0,0) and the other at (^1,^2). The behaviour shown in figure 2 agrees 
well with the qualitative discussion of the previous section. Some comments on the choice of 
parameters are in order. Unlike the case for large m^, negligible stop mixing [~ {Atip2—fJ''^i)] 
with two Higgs background fields necessarily implies negligible At and /x. Small fi may be 
problematic phenomenologically^. However, as we have mentioned, the transition takes place 
for low values of ttia along a field direction for which ip2/<fi is large. As a consequence, the 
effect of /i in the off-diagonal stop mixing is suppressed and its influence on the transition 
diminishes, which has been actually observed in [^. In practice, larger values of fi can 



be accommodated without weakening the strength of the transition. We can then safely 
set At = = for the purpose of analyzing the transition strength and this simplifies 
considerably the analysis. In the following we refer to tn and as mass eigenstates. 

Values of mq lower than 70 GeV would increase v/Tc, but they would be in conflict 
with Ap constraints |jT9|. Another possibility to increase v/T^ is to choose mfj < |]22|. 
In that case, tpt may get a vev and so has to be included in the discussion of the effective 
potential. Although, if realized, this region of parameters could have important implications 
(see [0, 0) we do not explore it in this paper and we stop at mu = 0. The perturbative 
analysis of the T 7^ potential along the squark direction (necessary to ensure that mfj < 
is allowed) may be problematic. We can make a rough estimate of the expansion parameter 
in the colour breaking minimum that may develop in the squark direction; in analogy with 
the SM, this expansion parameter would be the ratio of the squark quartic coupling over the 
SU{3) gauge coupling squared. This ratio is in principle of order one. In addition, when the 
thermally corrected mass of t/j is close to zero, radiative corrections can grow very large and 
affect the convergence of the perturbation series along the Higgs direction. These problems 
have been observed in ref. where a two- loop analysis of the potential (both along the 
Higgs and stop directions) was presented for the case mfj < 0. Clearly a non-perturbative 
analysis of this region {mfj < 0) would be desirable. 

However, sufficiently large values of v/Tc are obtained when two-loop corrections are 
taken into account [^, as we discuss in the next section, so that one is not confined to the 



choice mfj < 0. In particular, this means that the condition m^^ < nit is not required for a 
sufficiently strong phase transition. 



^We choose this too low value of Mt for the sake of the discussion as it corresponds roughly to the point 
were the effect of two-loop corrections on v/Tc is maximized. Results for Mt — 175 GeV are presented later. 
^We follow here the conventional definition of At and /i, which differs from the more general one used in 
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Figure 2: Order parameter v/T at Tc as a function of the pseudoscalar mass ruA- Curve (a) corresponds 
to the 2-loop resummed approximation, while (6) gives the 1-loop resummed result. Here Mt = 156 GeV, 
TTiQ — 70 GeV, mu — and tan/? — 2.5. 

5 Two-loop resummed potential 

For simplicity, we set g' = Q = hb 'm two-loop corrections. We only consider the case of 



negligible stop mixing, as this is the best case for a strong phase transition pO|, Q . 

The two-loop diagrams to consider are displayed in figure 3. Counterterm graphs are 
not shown. We represent Higgs bosons (A", if", , G^, G^, H^) by simple dashed lines, 
gauge bosons (g, W, Z, 7) by wiggly lines, squarks (tL,_R, by dashed lines with arrow, 
quarks (t,b) by continuous lines with arrow, and ghosts {cw, cz) by dotted lines with arrow. 
Within the stated approximations, the full two-loop corrections to the effective potential are 
collected in appendix B (we follow the notation of refs. 0). In the numerical analysis 



we use a high temperature expansion of the different integrals, keeping terms up to order 
gf, hf. In some corners of the parameter space we explore, for example the regions with 
large values of and/or ttlq, this high temperature expansion starts to fail. We control 
our numerical results in those cases interpolating between the region where the expansion is 
safe and the large m^,mQ region where the T 7^ contributions of the heavy Higgs and/or 
[i, b) L doublet are dropped. 

The dominant two- loop terms, i.e., those that have a greater influence on the transition 
strength, are written down in appendix C. These are terms involving logarithms of masses. 
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Cz) (z') (z") 



Figure 3: Two- loop vacuum graphs contributing to the effective potential (in Landau gauge). Dashed 
lines represent Higgs bosons, wiggled lines gauge bosons, dotted lines ghosts. Dashed lines with an arrow 
represent squarks, while quarks are solid lines with an arrow. 

In some cases, these logarithmic terms introduce a hnear dependence on the fields ipi and/or 
for small (pi^2 which is cancelled by other terms in the potential. We keep these terms 
to ensure that our approximation is well behaved in that respect. Finally, we also add the 
non-logarithmic terms which depend on the couplings Qs and ht and thus can be potentially 
large (although they mainly affect the transition temperature and do not change v/Tc too 
much) . 

The relative importance of the different diagrams follows the pattern discussed in ref. 
for the large case. For low values of the stop soft masses, the main effect is due to diagram 
(a') with gluon exchange. This QCD contribution significantly increases the value of v/T 
at the transition. This important enhancement is operative also in the low ttia regime as 
shown in fig. 2. Curve (a) gives v/Tc computed from the two-loop resummed potential as a 
function of the pseudoscalar mass m^. The corresponding one- loop result is given by curve 
(6). The parameters have been chosen to maximize the effect. 
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6 Results 



Before presenting the results of our numerical computation of the transition strength some 
comments on our definition of v/Tc are in order. We have taken as the temperature 
of degeneracy between the symmetric and broken-phase minima. A possible alternative is 
to choose the temperature Tq at which the symmetric minimum gets destabilized (that is, 
y'(0) = along some field direction); in that case, the necessary condition for baryogenesis 
is f /To Z 1.2 — 1.3. One can refine the precise number entering this condition by estimating 
the sphaleron energy [^, ^ . The true critical temperature should be in the narrow interval 
[To,Tc] if the transition is first-order. When two-loop corrections are included, in particular 
those of the form 5V ~ T^yj^logy), V"{ip = 0) diverges. This behaviour is associated 
with the (non-screened) transverse modes of gauge bosons. In principle, this prevents the 
computation of Tq. For practical purposes, one can overcome this difficulty by turning on a 
small magnetic mass (which, after all, is generated non-perturbatively^ and checking that 
Tq is not sensitive to its particular value. Although we have chosen to work with T^, which 
can be defined in a cleaner way, we checked that the results obtained for v/Tq are consistent 
with those that we present now. 



t.8 I 1 1 1 1 1 1 1 1 1 1 1 r 




0.8 I I ' l / I -^-^-1— r— , — I 

50 too 150 200 250 300 350 400 450 500 

rriA (GeV) 



Figure 4: Lines of v/Tc (solid) for different values of tan/? as a function of m^. The mass of the lightest 
CP-even Higgs is given by the dashed contour lines. Adt — 156 GeV , mjj = 0, ttiq = 70 GeV. 



^°In the abelian-Higgs model no such magnetic mass is generated, even non-perturbatively. However, the 
diagrams responsible for the divergence are also absent as was to be expected. 
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rriA (GeV) 



Figure 5: Same as fig. 4 but for Mt = 175 GeV, mu = 0, mg 250 GeV. 



In figure 4 we plot v/Tc (solid lines) as a function of the pseudoscalar mass tua for different 
values of tan (3. The other relevant parameters are Mf = 156 GeV (a low value that roughly 
corresponds to the maximal effect of two-loop corrections), m^/ = 0, and mq = 70 GeV 
(to satisfy the Ap constraint). In figure 5 we repeat the plot for Mf = 175 GeV {mjj = 0, 
rriQ = 250 GeV). These plots show that the transition is always stronger for lower values of 
tan P and larger values of m^. At large ttia, v/T^ gets stabilized and the results of ref. |^ are 
recoveredj^. At low tra, v/T^ drops and the curves for different values of tan/5 get focused 
in a narrow range. This behaviour, previously observed in refs. |26, 27], has been discussed 
in section 3. For a given value of tan/3, there is a critical value of itia below which the 
transition is too weak for baryogenesis. This critical value of itla increases with increasing 
tan/? and goes to infinity for tan/? ~ 3.5 — 4. For larger tan/5 the transition is too weak for 
any value of itla- 

Small tan/5 values (< 2) give a very strong transition because the Higgs mass irih is 
small. Dashed lines are contour lines for nth. In the way we are plotting our results, these 
curves are determined by the solid ones, because rrih depends on niA and tan (3. We use 
the LEP limits on rrih to determine the experimentally allowed region. For large m^, these 



^^The apparent discrepancy between figure 5 at large niA and the results of [ p3) for the same top mass 
is due to a more restrictive choice of mg in [ p3[ . There, a one-loop expression for Ap was used to estimate 
the lower limit on mq. This limit is weaker when QGD corrections are included (see |Q and references 
therein). 
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limits arise from the non-observation of the production process Z —>■ h^Z*, and we take the 
Standard Model lower limit rriho > 70.7 GeV, because in that regime has SM couplings. 
For light A^, the process Z — > h^A^ becomes important and in fact sets the strongest limit 
in the low m^-large tan/3 region. However, that region already corresponds to v/Tc < 1. 
For our purposes, we can take the line of mf^o = 70 GeV as the limit of the allowed region 
(it gets a few GeV weaker when tua ~ 100 GeV). After imposing the LEP bound we find 
from fig. 5 an absolute lower bound on rriA of about 120 GeV (below which the transition 
is too weak for baryogenesis)^ which corresponds to tan/5 ~ 3. This lower bound grows for 
larger tan/?. The upper limit on tan/? is tan/3 < 3.6 and the lower limit is tan/3 > 2.25. The 
Higgs mass is below 85 GeV, within the reach of LEP II. We find no trace of a good region 
for baryogenesis at low rriA and large tan (3. 

Our results are in good agreement with those obtained in ref. using a 3d reduced 
effective theory and lattice Monte Carlo results. The upper bound on rrih of about 70 GeV 
quoted in that paper is smaller than our bound of 85 GeV because of the choice = 
100 GeV in [^. Smaller values of mu would allow for larger Higgs masses, as is shown in 
figure 5 of that paper, but in that region the approximations used to construct the reduced 
3d theory start to break down. The same comments apply to the lower bound on niA of 
tua ~ 200 GeV, which can be lowered for smaller values of mu- 

An independent 3d analysis of the transition has been carried out in refs. ETf. It is more 



difficult to compare our results to theirs because they imposed SUGRA-type constraints on 
the parameters of the model. Of course, this restriction leads to stronger bounds on ruh, tua 
and tan /3 but the qualitative dependence of the transition strength on different parameters 
is similar to our results. 



7 Conclusion 

In this paper, we have studied the electroweak phase transition in the MSSM, searching for 
regions where it is strong enough to permit electroweak baryogenesis. If stops are only mod- 
erately heavy they control the strength of the transition and the use of perturbation theory 
for its study is well justified. We first addressed the question of whether the electroweak 
phase transition can be a crossover in the MSSM. We clarified the differences with respect to 
the Standard Model case and estimated the critical Higgs mass for crossover, which depends 
on how stops affect the transition. This critical Higgs mass (beyond which the analytic 
crossover can take place) can be significantly larger than in the Standard Model. 

We scanned the parameter space {m a, tan (3), qualitatively explaining the infiuence on 
the transition strength expected from different effects. We identified a region with light 
(unmixed) stops, light rrih, and small tan/? where the transition can be strong. That region is 

^^For low values of tua, the charged Higgs boson is also light and gives a large contribution to 6 ^ 57 
which may conflict with the measured value of this branching ratio. However, cancellations with other 
supersymmetric contributions (from light stops and charginos) may occur and we cannot exclude these low 
values of without further assumptions on the SUSY parameters. 
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widest for large m^, and shrinks when ttla is lowered. We then focused on the interesting low 
irtA region and performed a numerical study evaluating the T-dependent effective potential 
of the model (simplified to a two Higgs doublet model plus third generation squarks) up to 
two-loop resummed contributions. We found that the region where the transition is strong 
extends down to ~ 120 GeV with roughly 2 < tan /3 < 4 and nih ~ 85 GeV (light stops 
are required but mi < rrit is not a necessary condition). This region will be explored by LEP 
II in the near future, which will stringently test the viability of electroweak baryogenesis in 
the MSSM. 
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A Field-dependent Masses and Mixing Angles 

We give here the field-dependent masses (in the background (pi^2 = (-f^i2)v^) of different 
species of particles relevant for the effective potential. Where appropriate, field-dependent 
mixing angles are defined. We also give the leading parts of thermal masses, needed for 
resummation of the potential. These masses are computed assuming that the plasma is 
populated by SM particles, an extra Higgs doublet and squarks of the third generation. 

• Gauge bosons: 

The number of degrees of freedom is n^y = 6 (with ti-Wl — '^i'^Wt = 4) and nz^^ = 3 
{nzi^,-yj^ = l,nzj.,-yj. = 2). The masses are 



M^ = ^,V + ^^), M^3_^ = i(_{, + (A.1) 



which gives 

Ml = 

4 



MI = Ig'{^1 + ^1), M' = 0, (A.2) 



where = g'^ + g'"^. 

At finite temperature the leading contribution to the thermal mass is zero (at leading 
order gT and to all orders in perturbation theory) for transverse modes while for 
longitudinal modes the above expressions change to 



M(V3-B). - M^3_^ + ( ^ ) , (A.3) 



with Ub^ = y|(7'^T^. The eigenvalues of M^^^^_^^^ are and M^^. In the approxi- 
mation g' = (used for the two-loop corrections in this paper) it follows that M^^ = 
and Mz^ = Mwr^ = Ml- 

For gluons {ug = 3x8) the masses are 



Ml=0, M'^^=0 + UgL = lg'sT'. (A.4) 



3^ 



Quarks: 

For the third generation we have {rit = nf, = — 12) 



rut = -^ht(f2, rrib = -^hbcpi. (A.5) 



where ht {ht) is the top (bottom) Yukawa coupling. 
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Squarks: 

We consider only the third generation, that is, stops and sbottoms (with n^^ — n^^ — 
= — 6) . The mass eigenstates are given by the diagonahzation of their 2x2 
mass matrix defined in terms of interaction eigenstates. For stops this is given by: 

V ^ht{Av2 + l^Vi) ml + ml + \g'\^l - ^1) )' 

while for sbottoms it is: 

I ml + ml-l^{3g^ + g'^){^i-^l) j^h{Ai,Vi + fi^2) \ 
\ -^h{Ab^i + fi^2) rnl + ml- ^g'^{^l-^l) ) - 

At finite T, the masses corrected by thermal effects are given (at leading order) by 





Ml = Ml+\ - ] , (A.6) 



with 



2t-t2 I 4 /2t-t2 I -'- i,2t-t2 



g^S 3 

= IglT' + Yl'"^' + \hlT\ 

Mixing angles can be defined in the usual way. As we are interested in the case of 
negligible squark mixing we do not give expressions for these angles explicitly. In this 
case the mass and interaction eigenstates coincide. 

Higgs sector: 

The two Higgs doublets are 



V2 



The background </7i,2 is CP conserving so that the scalar {h^,h!^), pseudoscalar 
(/i?\ h!^) and charged (-fff*, H2) sectors can be treated separately. For each of them 
we have a 2 x 2 hermitian mass matrix 



,,2 / mw;^ 1 



(where k — r,i,c) and a mixing angle dk which relates mass eigenstates {H,h) to 
interaction eigenstates {Hi, H2) according to 



iff) j \ sin^fe cos^fc 




Hi ' \ _ I cos9k - sin 6'^ , , , 
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The mixing angles are then given by: 
sin 26k = 



12 



(A.IO) 



cos 26k 
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For the CP-even sector, ifj^ = ^?2, ^^''^ = ^^''^ = "-/.o = 
reads: 

Z!r.W_,o2^ rr^2 G2 

and, at finite T, 



1 and 



l(r) 



with 



4^ 12^ 



i^2y2 ^ 1/2^2 ^ 3 2 2 

4^ 12^ 4 * 



(A.ll) 



(A.12) 



(A.13) 



There is a mixing angle 6r for the diagonalization of ([A.ll|) and a T-dependent angle 
6r for the diagonalization of (|A.12| ). 

For the CP-odd sector, H^l = h^]^, h^"^ = G°, H^^ = Ugo = n^o = 1 and (Q 
reads: 



— m; 



8 VV-l r2y '"'12 



(A.14) 



and, at finite T, an equation similar to ( A.12 ) holds, with r —>■ i. 

For the charged sector, hI"'^ = H{*, H^^ = , h^^^ = G+, = H+, nG± = nH± 
2 and ( [Al8|) reads: 



8 'f'l + 8 ^2 



-"^?2 + T'^1'^2 



(A.15) 



^1 ' +^-^ 

V -"^?2 + T'^1'^2 

and, at finite T, an equation similar to ( |A.12| ) holds, with r ^ c. 

For different cross-checking purposes it is useful to consider the behaviour of the mixing 
angles 6k in two different limits: 

(a)\n the large limit 



6r 



TT 



6i, 6c 



TT 



(A.16) 



(b) In the physical limit, ^91,2 — > ^^1,2, 

6r ^ a 



01, dc 



(A.17) 
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where a is the physical mixing angle in the CP-even sector. 

For later use it is convenient to define the following abbreviations (with the indicated 
limiting values [(a), (6)]): 

{f>Q = ipi cos 6r + (p2 sin 9r 
= COS 9r — (pi sin 9r 

(fsin = (fl COS 6r - (p2 siu 6r 

{pcos = V^2 COS 9r + (pi sin 9r 
(fQ^i = (p2 sin 6i - (fi COS 6i 
(fy^i = (p2 COS 6i + (fi sin 6i 
V^cos.c = f2 COS 9c - ipi sin 9^ 

^sin,c = ^2 sin 9c + ipi COS 9c 



[0, V cos(p — a)\ 
[(p, V sin(/5 — a)] 
[ip sin 2/3, V cos(/5 + a)] 
[(p cos 2/3, t> sin(/3 + a)] 
[0,0] 

[— V9COS 2(3, —vcos'. 
[y3sin2/?, f sin 2/3]. 



(A.18) 



B Two-loop Resummed Potential 

We present the resummed two-loop corrections to the finite temperature effective potential 
in the MSSM framework described in the text. We make the approximation g' = hfy = 0, so 
that in the formulae below we make no distinction between mz and mw which are called M 
[Ml is the longitudinal component). We further assume that left-right mixing in the squark 



sector is small and can be neglected. We follow the notation of p8| , ^ and the labeling 
of different contributions corresponds to our figure 3. As usual, Nc = 3 counts the number 
of colours and we use the abbreviations Cr = cos 9r, S2c = sin2^c, etc. For some recurrent 
combinations of fields and angles (like e.g. ipiCr + ^2Sr) we use the short-hand expressions 
defined at the end of appendix A. With our resummation method in the scalar sector (we 
included the thermal mass correction in all Matsubara modes) masses and mixing angles are 
T-dependent: mj, 9k- For simplicity we drop the bars in the formulae below and simply 
write mj, 9k everywhere. 

- ^—[2cos\9, + 9c)[Vssv{mG±,mh»,M)+Vssv{mH±,mHO,M)] 

+ 2 sin^ {9r + 9c) [T>ssv ("^g± , , M) + Vssv ("^/f± , , M)] 

+ 2 cos^ {9c - 9i)[DssvimG±,mG0, M) + Vssv{mH±,mAo, M)] 

+ 2sin^(6'c - 9i)[VssvimG±,mA0,M) + VssvimH±,mGO, M)] 

+ cos'^{9i + 9r)['Dssv{'mGO,mhO,M) + DssvimAo^niHO, M)] 

+ sm^{9i + 9r)[Vssv{mG0,mH0,M) + Vssv{rnhQ,mAo, M)] 

+ Vssv imG± , mG± , M) + Vssv imH± , , M)] } , 

0"')= - ^Nc[Vssv{mi^,mi^,M)+Vssv{rn-,^,rn-,^,M)+4Vssv{rni^,rn-,^,M) 



- ^{N^c - 1) [Vssvim^,,mi^,0)+Vssvim-,^,m-,^,0) 



.2 



4 
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— 

yip) ^ 



- — / [(^^P5yy(mfe0, M, M) + iplVsvvirriHO, M, M)] 

- [2Vsv{mG±,M) + Vsvimco, M) + D^yl^M', M) 



AglVffvimt, rut, 0) - K, m^, M) - P//y(0, 0, M)] 

0, M) - Vffv{Q, 0, M)] - ?,g^ntVf^v% 0, M), 

4 * 



clVffs{mt, mt, rriho) + s^Vffsimt, rrit, mno) + cfVffs{mt, rrit, mao) 



+ s1Vffs{mt, mt, rriAo) + 2clVffs{mt, nib, mG±) + 2slVffs{mt, nib, mH±] 

- [Dvvv{M, M, M) + ?,Vllt{Ml, Ml, M) - Wllt{M, M, M)] , 

- 3g''Vr„v(M), 

- -^g^Vyy{M,M), 

- -^g^ {2 - S2c</?5m)^H(mfto, mG±,mG±) + {(fv + S2c<fsinfiiimho,mH±,mH±) 



+ 2 



[i^o - S2c</7co5)^H(m//o, m^i, m^i) + {^Po + S2cVcosfY{{mH'^,mH±, mH±) 



+ 4c: 
+ 4 
+ 3c; 



2c 



ipliY{{niGO,niG±,niH±) + iB.{niAo,mG±,mH±) 



2 TT/ 



+ (2(/?o - 3c2r</?sm)^H(mj/o, m^o, m^o) + (2(/7^ - 3c2r</?cos)^H(mijo, m^fo, m^o) 

+ ^'/'L [H(m^o, mco, mco) + H(mhO, 771^0,771,^0) 

+ cli(flin [H (777^^0, 777^0, mco) + H(777i70 , 777^0 , 777 AO ) 

+ 24 [v?LH(mfeO, 777 G0, 777^0) +95LH(mHO, 777GO,mAo)]}, 

y^') = - ^Nc{[htCrip2 - t/Vcos/4]^H(777ft0 , 777^-^, 777^^) + {h^Cr(p2fli{mhO , m^^, 777^^) 

+ [/7tSr-</?2 + 5'Vsm/4]^H(777ifO, mi^,mi^) + {htSrip2fB.{mH0,mf^,mfJ 

+ [htCc(p2 + g'^(psin,c/'^?ti{mG±:mt^,miJ + [g'^(pcos/^?ti{mho, mi^,mi^) 

+ [hlSc(p2 - fi'Vcos,c/2]^H(777^/±, 777^^, 777^^) + [fif V^im/4]^H(777ijO , 777^^ , 777^^ ) } , 
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16 

+ 2 [(1 + S2cS2r) [Vs{mG±,mho) + Vs{mH±,mHo)] 

+ (1 - S2cS2r) [VsimG±,mHo) + V s(mH± ,mho)] 

+ (1 - S2cS2i) [Vs{mG±,mG») +Vs{mA»,mH±)] 

(1 + S2cS2i) [Vs{mG±,mAo) +'Ds{mGO,mH±)]] 
3 



+ 
+ 



cl^ ['Ds{mho,mho) + Vs{mHO,mHo)] + cl^ ['Ds{mGo,mGo) + 'Ds{mAo,mAo)] 
+ (2 - 3cl^)Vs{mhO,mHo) + (2 - 3cli)Vs{mGO,mAo) 

+ C2i C2r [Ds {rriho , mco ) + Vs {niHO , m ao )-T>s{mho, m ao ) - Vg {niHO , mco ) ] } , 
+ [T^s {mtL ' "^G° ) + '^s {rni^ , mco )] + s1[Vs {mi^ , m ) + Vs {mi^ , m^o ) ] 

+ 2c^[^5(^?«,^G±) +^5(mSi'"^G±)] +2s^[p5(^t«,^H±) +^5(^Si'"^i?±)]} 

+ ^fl'^-^c {c2r[C>s{mi^,rnho) - Vsimi^,mho) + X>s(mf^, m^fo) - Vs{mi^,mHo)] 

+ C2i [Vs {mi^ , ) - Vs {mi^ , 771^0 ) + Vs (m^^ , men ) - "P^ ("^t^ , "^g" )] 

+ '^C2c[Vs{mi^,mH±) -Vs{mi^,mH±) +Vs{mi^,mG±) - Vs{mi^,mG±)\ \ , 

iVe(7Vc + 1) m,-J + Vs{m-,^, m-,^ 



+ 
+ 
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N,{N, + 1) Vs{mi^,miJ +Vs{m-,^,m-J 



(B.l) 



where 



Vs{mi,m2) = /(mi)/(m2). 

In the previous formulae, dimensional regularization (with n — 1 = 3 — 2e) is used to 
evaluate divergent integrals. Poles in 1/e and t^-dependent terms cancel when counterterms 
are included. In addition, counterterms contribute a finite piece to the potential. 

With g' — the T — counterterm potential reads 



SVr 



1 r3 



count — 



2 U 



■(3 - 2e)/^ {6Z,2 + SZ^2 + SZ^.) g\^l + ^l) - 6h'j}f,{6Z,2 + 6Z^2)ifl\ , 



with 



rp2 

l0 = J^{^ + ^^e 



■^[l + e(.,-log4)]. 



The MS renormalization functions 5Z (calculated in a model with two Higgs doublets 
plus third generation squarks) are 



A7 ^ 2 1 



24 



^Zh2 = 

The finite contribution is then 



SVr 



count, fin 



64712 



M^Sh'^t - 2g^) + m^t{16gl - Zh]) log 2 



(B.2) 



The resummation procedure we followed for scalars can be implemented by adding and 
subtracting the appropriate thermal mass terms for the scalars in the Lagrangian. In this way 
the unperturbed Lagrangian contains already thermally corrected masses and thermal coun- 
terterms appear. To cancel all divergences, diagrams involving those counterterms should 
be included. For the two-loop potential, one-loop thermal counterterm contributions are 
needed and give 

5Vth,count = "2 ^inj7(m-), (B.3) 

i,scalars 

where Ilj is the thermal self-energy of the i*'^ scalar. 



C High T expansion of the Potential 

We give the dominant two-loop terms of the potential presented in appendix B using a high 
T expansion. 

First we have logarithmic terms coming from bosonic setting-sun diagrams: 

^2 ( 



^ T^<!cos2( 



647r2 



Oc) 



+ sin^(6'r + 9c 



+ cos" {9c -9i) 



+ 



+ 



+ 



+ 



(M^ - 2Tr4± - 2m^o)log 



3T 



mG± + rriifi + M 

mG± + ruho 



M2 



log ( 

\mG± + ruho + M 



+ (M^ - 2m|± - 2m^o) log 



3T 



m 



2 \2 



M2 



log 



(M^ - 2mG± - 2m]jo) log 



3T 



mG± + TTlHO + M 



m 



2 ^2 



log 



/ m(3± + rriHO 



3T 



+ (M^ -2m^± -2m^o)log 



M2 



■log 



{M^ - 2m%± - 2mGo) log 



3T 



mQ± + mQO + M , 
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+^{M^ -Aml±)log 



3T 



2mG± + M 



log 

+2 log 
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m,,() + 2M 
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+ 
+ 
+ 



+ 
+ 
+ 



+ 
+ 
+ 



2M4 
+ Slog 



+ 
+ 
+ 
+ 

log 



{mQ± - mQo) / mG± + moo 
M2 \mG± + moo + M 

f 3T 

-2m^± -2m^„)log 



■log 



M2 

- 2mG± - 2m^o)log 



3T 



(C.l) 



m^i + m^o + M 



- 2m%± - 2mro) log , 

Vmjf± + moo + M 



M2 

-2m^o -2m^o)log 



\mif± + 777.(50 + M 

3T 



M2 



log 



mfjo + rriho 



- 2m\o - 2mjjo) log 
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■log 
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3T 



i^ojiOg 
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M2 
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-2m\o -2777^0) log 
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3r 



-(M^ - 4777|^±)log 
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j +log 
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+2 log 



/ mffo + 2M 
\mH0+2ML 



+ Slog 



3T 



ruffo + 2M 



+ 



M2 



log 
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33M2 log — - 6(M2 - 4M^) log , 



KniHO + 2M. 
3T 



(C.3) 
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3T 



rufio + 2ma± 



+ 2(930 - S2c^cosf log 
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3T 



+ 2{(po + S2cVcosy log 
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3T 



rriHo + 2mH± 



+ 2((/7^ + S2c</'sm)^log 

\m/iO + 2mH± 
2 , 3T \ 2 ^ ( 3r 

V'sin log + ^cos log 

3T ^ 2 ■ ^ 3T 



+ 4c; 



2c 



+ 4 



<^o,ilog 



moo + m(3± + mH± 
3T 



V'^.i log 



+ (2(/7o - 3c2,.(/?sin)^ log 



+ (2(/7^ - 3c2,.(/?cos)^ log 



3T 



log 



3T 



log 



+ 'Psin 

+ Hi 

and including squarks 
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+ log 



+ log 



3T 



. 1 - 3T 
'/'co.log 



3r 



m/jO + moo + m^o 



+ </'Llog 
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^ T^AT, <! 4(M2 - 2m? - 2m? ) lo; 
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mi^ + m^^ + 
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(m? — m? Y 
^log 



3r 



+ (M^ - 4my log [ ^^.^^m ) + " 



3r 



2m^^ + 



m^-i)|<iog(;^)+<i°g 



3T 



167r2 



+ log ' » i ™2 



.2mf^> 

3T \ 21 / 3T 
+ "^6^ log 
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2m 
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N, 
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TU{htCrV^2r log 



3T ' 
mho + 2mi^^ 



+ {htSr(p2flog 



3T 
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+ {htCrCp2 - -r^cos) log — + {htSrip2 + -T^sin) log — 

4 \mhO + 2mt^J 4 \mH0 + 2mf^J 

+ {^^cosf log ( ) + {^^sinf log ( ) 

+ {KCc^2- ^^cos.c) 

2 \mG±+mi^+miJ 

+ {hfSc^2 - ^^sin,cf log ( . . J I (C.6) 

The linear terms needed to insure dV/d(pi = dV jd^p^ = at = = are 
3 

^^lin = --^9'^M{mho + 2mG± + moo + uiho + 2mH± + "^^o)t^> (C-T) 



(C.8) 



(C.9) 



^2 



167r2' 



and 



3M 



j'2 



3g^NcM{mt^ + m, 



T2 



It can be easily checked that (i = 1, 2) 

d 



dep. 



Hn 



327r2' 



= 0, 



¥3=0 



a similar relation for (a' + c') , and 

d 



(6+n) 
log 



0. 



The non-logarithmic terms involving only Qg and /ij couplings are: 
For Standard Model diagrams 



647r2 



16^2 ( ^ log 2 - 1 - cb) + 37Vc/i? log 2 - cs^ 



(C.IO) 

(C.ll) 
(C.12) 



(C.13) 
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with Cb = log(47r) — 7e, and from (|B.2| 



5V,ount = (16^,2 - 3/1^) log 2. 



647r2 



(C.14) 



Diagrams involving squarks give 



647r2 



(iV,^-l)(c2-l)(m4+m,^J, 



327r2 



(iV2 



n 



C2hf(pl 



1287r2 



167r2 



6 



+ 5,,V^^"\ (C.15) 



with C2 ~ 3.3025. 

The terms dcgV add up to 



i, scalar s 



The contribution from thermal counterterms, eq. (p.3| ), gives exactly the same but with 
opposite sign so that they cancel out. If resummation of scalars is implemented only on zero 
Matsubara modes, there are no thermal counterterms. In that case 5cgV combines with the 
one-loop unresummed scalar contribution 

Cb 



to give 



327r' 



Cb 
327r^ 



(C.16) 



(C.17) 



That is precisely the one-loop result if we resum all modes, so that after the cancellation 
of two-loop contributions and thermal counterterms, both resummation methods give the 
same result. This is no longer true in the presence of particle mixing (this occurs in our case 
in the Higgs sector) where a small numerical difference is expected. 
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